The physical parameters, excitation and chemistry of the rim, jets 
and knots of the planetary nebula NGC 7009 ^ 



Denise R. Gongalves 
Instituto de Astrofisica de Canarias, 
^f) • E-38205 La Laguna, Tenerife, Spain 

' e-mail: denise@ll.iac.es 

o ■ 

Q_i! Romano L. M. Corradi 

O ' Isaac Newton Group of Telescopes, 

Apartado de Correos 321, E-38700, Sta. Cruz de la Palma, Spain 



00 



0^ ' e-mail: rcorradi@ing.iac.es 

^ I Antonio Mampaso 

■ Instituto de Astrofisica de Canarias, 

I E-38205 La Laguna, Tenerife, Spain 

^SJ ■ e-mail: amr@ll.iac.es 

o _ 

^ ■ Mario Perinotto 

, Dipartimento di Astronomia e Scienza dello Spazio, Universitd di Firenze, 

• Largo E. Fermi 5, 50125 Firenze, Italy 

e-mail: mariop@arcetri.astro.it 

o 



ABSTRACT 

We present long-slit optical spectra along the major axis of the planetary nebula NGC 7009. 
■ These data allow us to discuss the physical, excitation and chemical properties of all the morpho- 

I logical components of the nebula, including its remarkable systems of knots and jets. The main 

results of this analysis are the following: i) the electron temperature throughout the nebula is 
remarkably constant, Te[0 in] = 10 200 K; ii) the bright inner rim and inner pair of knots have 
similar densities of ~ 6000 cm~'^, whereas a much lower density of ~ 1500 cm~^ is derived 
for the outer knots as well as for the jets; Hi) all the regions (rim, inner knots, jets and outer 
knots) are mainly radiatively excited; and iv) there are no clear abundance changes across the 
nebula for He, O, Ne, or S. There is a marginal evidence for an overabundance of nitrogen in the 
outer knots (ansae), but the inner ones (caps) and the rim have similar N/H values that are at 
variance with previous results. Our data are compared to the predictions of theoretical models, 
from which we conclude that the knots at the head of the jets are not matter accumulated during 
the jet expansion through the circumstellar medium, neither can their origin be explained by the 
proposed HD or MHD interacting-wind models for the formation of jets/ansae, since the densities 
as well as the main excitation mechanisms of the knots, disagree with model predictions. 

Subject headings: planetary nebulae: individual (NGC 7009) - ISM: kinematics and dynamics - ISM: jets 
and outflows 
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1. Introduction 

In addition to their large-scale structures, 
mainly identified in the forbidden [OIII] emis- 
sion line, many planetary nebulae (PNe) have a 
number of smaller-scale structures which are in- 
stead more prominent in low-ionization lines, such 
as [N ii], [O li] and [S ii]. In a previous paper 
(Goncalves, Corradi, & Mampaso 2001), we have 
used the acronymous LISs to identify these low- 
ionization structures. LISs appear with different 
morphological and kinematical properties in the 
form of pairs of knots, filaments, jets, or isolated 
features moving with supersonic velocities through 
the large-scale components in which they are lo- 
cated, or instead as structures with the above 
morphologies but with low velocities that do not 
differ substantially from that of the main nebula^. 
LISs can be easily seen in the imaging catalogues 
of Balick (1987), Schwarz, Corradi, & Melnick 
(1992), Manchado et al. (1996), Corradi et al. 
(1996) and Corny et al. (1999). For a detailed 
study concerning the morphological and kinemat- 
ical classification of LISs in PNe, see Gongalves et 
al. (2001). 

A very interesting class of LISs are jets. These 
are highly coUimated outflows (with aspect ratios 
varing between 3 and 20; Lopez 1997), which are 
directed radially outward from the central star, 
appear in opposite symmetrical pairs, and move 
with velocities substantially larger than those of 
the ambient gas that forms the main bodies of the 
nebulae. It should be noted in addition that jets 
are much more collimated than other bipolar out- 
flows of PNe, i.e, more collimated than the lobes 
of a bipolar PN. 

NGC 7009, the "Saturn Nebula", is a well stud- 
ied elliptical PN, which possesses a jet-like sys- 
tem as well as two pairs of low-ionization knots 
along its major axis. A number of previous stud- 
ies showed that high-excitation lines dominate the 
inner regions along the minor axis, while the low- 

^Bascd on observations obtained at the 2.5 m Isaac New- 
ton Telescope (INT) of the European Northern Observa- 
tory, and with the NASA/ESA Hubble Space Telescope, 
obtained at tlie Space Telescope Science Institute, which is 
operated by AURA for NASA under contract NAS5-26555. 
^The fast, low-ionization emission regions, FLIERs (Balick 
et al. 1993) and bipolar, rotating, episodic jets, BRETs 
(Lopez, Vazquez, & Rodriguez 1995) are particular types 
of LISs, i.e., high- velocity knots and jets. 



excitation ones are enhanced at the ends of the 
major axis. The ionization structure is further 
enriched by the fact that the low-excitation re- 
gion present strong variations in excitation level 
and dumpiness. NCC 7009 was classified as an 
oxygen-rich PN (Hyung & Aller 1995), with an 
0/C ratio exceeding 1, and anomalous N, O, and 
C abundances (Baker 1983; Balick et al. 1994, 
Hyung & Aller 1995). Its central star is an H-rich 
0-typc star, with effective temperature of 82 000 K 
(Mendez, Kudritzki, & Herrero 1992; Kingsburgh 
& Barlow 1992). The kinematics of NGC 7009 
was studied by Reay & Atherton (1985) and Bal- 
ick, Preston & Icke (1987). The distance deter- 
minations of NGC 7009 vary from 0.5 to 2.3 kpc 
from various indirect methods (Acker et al. 1992; 
Cahn, Kalcr & Stanghellini 1992; Maciel 1995). 

According to the theoretical models for the for- 
mation of jets and of the pairs of knots that are of- 
ten found at their tips (Soker 1990; Garcfa-Segura 
1997; Garcfa-Segura et al. 1999; Garcia-Segura & 
Lopez 2000; Steffen, Lopez, & Lim 2001: Gardiner 
& Frank 2001; Blackman, Frank, & Welch 2001), 
a significant density contrast between the jet it- 
self and the knots at its tip is expected. This is 
a quantity that can in principle be determined by 
optical spectroscopy. 

In addition, previous observational work sug- 
gests that the LISs of NGC 7009 and other PNe 
have higher abundances of nitrogen than their 
main bodies (see Balick et al. 1994 and references 
therein). These authors proposed that such an 
overabundance might indicate that the knots of 
NGC 7009 were ejected at a different epoch than 
the main body of the nebula, when the stellar sur- 
face was enriched of these elements due to the 
drastic changes of nitrogen in AGB stars and some 
dredge-up episode, or that the knots come from 
chemically peculiar regions in the stellar surface 
or its interior. 

In this paper, we test the above two issues by 
obtaining spectroscopic data for NGC 7009, de- 
riving its physical and chemical properties, and 
finally showing that neither a density contrast be- 
tween the jets and knots is found, nor are their 
abundances strongly enhanced compared to the 
rest of the nebula. In the following sections, we 
first describe the data, their reduction and the var- 
ious morphological/kinematical structures of the 
nebulae (Sections 2 and 3). Then we proceed with 
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the analysis of the PN spectra in terms of phys- 
ical, excitation and chemical parameters (Section 
4) and, finally, we discuss the implication of the 
present results for the formation of the LISs of 
NGC 7009 and give our conclusions, in Sections 5 
and 6. 

2. Observations and Data Reduction 

Spectra of NGC 7009 (PN G037.7-34.5) were 
obtained on 2001 August 29 at the 2.5 m Isaac 
Newton Telescope (INT) at the Observatorio del 
Roque de los Muchachos (European Northern Ob- 
servatory, La Palma, Spain), using the Intermedi- 
ate Dispersion Spectrograph (IDS). The IDS was 
used with the 235 mm camera and the R300V grat- 
ing, providing a spectral coverage from 3650 A to 
7000 A with a spectral reciprocal dispersion of 3.3 
A pix~^. The spatial scale of the instrument was 
0".70 pix-i with the TEK5 CCD. Seeing varied 
from 0".9 to The slit width and length were 

1.5" and 4', respectively. The slit was positioned 
through the center of the nebula at P. A. = 79°, 
and the exposure times were 3x30 s, 3 x 120 s, 
and 3 x 600 s. 

During the night, bias frames, twilight and 
tungsten flat-field exposures, wavelength calibra- 
tions and exposures of standard stars (BD-l-332642, 
Cyg OB2#9, HD19445, and BD-h254655) were 
obtained. Spectra were reduced following the 
IRAF instructions for long-slit spectra, being 
bias-subtracted, flat-fielded, combined in order 
to improve the signal-to-noise ratio and elimi- 
nate cosmic rays, wavelength calibrated, and sky- 
subtracted. Finally, they were flux-calibrated us- 
ing the above-mentioned standard stars and the 
mean atmospheric extinction curve for La Palma. 

We also retrieved HST [N ii] and [O III] archive 
images of NGC 7009 obtained on 1996 April 
28, and on 2000 April 7, respectively, with the 
WFPC2 camera in the F658N and F502N filters, 
with total exposure times of 1200 s and 320 s. 

3. The morphological and kinematical 
components of NGC 7009: previous 
v^rorks 

The HST [O III] and [N ii] images of NGC 7009 
are shown in Figure 1. The nebula consists of sev- 
eral components. The brightest one is the ellipti- 
cal inner shell {rim) roughly elongated in the EW 




Fig. 1.— The HST [O ill] and [N ii] images of 
NGC 7009 on a logarithmic intensity scale. Slit 
position is indicated by short lines (P. A. 79°), 
while labels mark the position of the several struc- 
tures, where "K" , "J" and "R" stand for "knot" , 
"jetlike" , and "rim" , respectively. Extensions of 
the selected structures with respect to the central 
star are: from -26.6" to -21.7", Ki; from -20.3" 
to -16.1", Ji; from -14.7" to -6.3", K2; from 
-4.9" to -1.4", Ri; from 2.1" to 3.5"; R2; from 
4.9" to 14.0", Kg; from 17.5" to 22.4", J2; and 
from 23.8" to 28.7", K4. 

direction, where it extends out to ±13" from the 
center. The inner portions of the rim cut by our 
spectroscopic slit are indicated in Figure 1 (bot- 
tom) by the labels Ri and R2. 

Surrounding the rim, a fainter elliptical outer 
shell {shell) that shares its main axis with the rim 
is also visible in the HST images. In [O ill], at 
fainter intensity levels the shell shows two "polar" 
bubbles. Inside these bubbles, coUimated LISs are 
clearly seen in the [N li] images. These are com- 
posed of a "pencil" of emission (the jetlike struc- 
tures, labeled Ji and J2) that connects the inner 
rim with two outer knots (Ki and K4). The [N 11] 
image also shows another pair of microstructures, 
which appear as "compact groups of bright knots" 
that we call the inner knots (K2 and K3). The ex- 
tension of all features is given in the caption of 
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Figure 1. 

The rim of NGC 7009 is likely to be formed 
by the interaction between the fast post-AGB and 
the slow AGB winds (Kwok, Purton, & Fitzgerald 
1978; Frank & Mellema 1994; Mellema & Frank 
1995; Zhang & Kwok 1998, Schonberner 2002). 
The cavity delimited by the rim is filled with a 
hot (T = 1.8 X 10*' K) gas emitting thermal X-rays 
and its size is in agreement with that of the optical 
rim (Guerrero, Grucndl, & Chu, 2002) providing 
a direct proof for the shocked fast stellar wind in 
the interacting winds scenario. The formation of 
the shell has been ascribed to the action of the 
photoionization front on the AGB matter not yet 
reached by the shock produced by the fast wind 
(e.g., Schonberner 2002). 

The origin and nature of the LISs in NGC 7009 
are rather unclear at present. The inner and outer 
pairs of knots have been studied by many authors 
(Reay & Atherton 1985; BaHck et al. 1994; Hyung 
& Aller 1995; Lame & Pogge 1996; Balick et al. 
1998), who often called them caps and ansae, re- 
spectively. Reay & Atherton (1985) have stud- 
ied the kinematics of the inner and outer knots 
m the [O i] 6300 A line using Fabry-Perot data. 
They derived an expansion velocity for the inner 
knots of ±38kms~^ with an inclination with re- 
spect to the line of sight i = 51°, and a velocity of 
±60 kms""'^ with i = 84° for the outer pair. Note 
that the determination of the deprojected veloc- 
ities and inclination angles by Reay & Atherton 
(1985) is mainly based on the proper motion mea- 
surements of Filler (1965), as the velocity in the 
plane of the sky appears one order of magnitude 
larger than the Doppler shift. Wo have estimated 
the separation of the outer knots in the [N ii] HST 
image taken in 1996, finding ^50.8 ± 0.5 arcsec, 
the uncertainty depending on the highly irregular 
morphology of these knots that makes it difficult 
to define a centroid. This is clearly not consistent 
with the extrapolation of Filler's measurements 
(separation of 48 arcsec in 1965 and apparent ex- 
pansion of 1.6 arcsec per century for each knot), 
making the velocity and inclination figures of Reay 
& Atherton (1985) quite uncertain. Their analy- 
sis also depends on the distance adopted, which is 
poorly known, as mentioned in the introduction. 

No velocity information is available from the 
literature for the jetlike structures connecting the 
rim to the outer knots. If we assume a continu- 



ity between the velocity of the outer knots and 
these jetlike structures, then the latter might be 
real supersonic jets. A more detailed study of the 
kinematics of NGC 7009, and in particular of the 
faint "pencil" and large-scale structures surround- 
ing it (shell and polar bubbles), is clearly needed 
to understand its nature better. 

4. Data Analysis 

Line intensities were mainly measured from the 
deep spectrum with 3x600 s exposure time. The 
shorter-exposure spectra were used to measure the 
bright H/3, [O iii]. Ha, and [N ii] lines that are 
saturated in the deep spectra. Measurements are 
given in Table 1, which lists the wavelength of the 
line identification and the corresponding ion (col- 
umn 1), the line fluxes measured for the different 
spatial regions selected (columns 2 to 9, as defined 
in Fig. 1 and Section 3), and the flux integrated 
along the slit for the whole nebula, labeled as NEB 
(column 10). All fluxes are normalized to /(H/3) 
= 100 in each of the regions considered. 

Errors in the fluxes were calculated taking into 
account the statistical error in the measurement 
of the fluxes, as well as systematic errors of the 
flux calibrations, background determination, and 
sky subtraction. Table 2 gives the estimated ac- 
curacy of the measured fluxes for a range of line 
fluxes (relative to H/3) in each of the selected re- 
gions. Absolute H/3 fluxes integrated along the slit 
in each region are also given in Table 2. 

4.1. Extinction 

Fluxes were extinction-corrected by using 
ch/3=0.16± 0.01 (the logarithmic ratio between 
observed and dereddened H/3 fluxes), determined 
from the Ha/H/3 ratio in the brightest regions of 
the nebula (Ri and R2). For the derivation of 
ch/3, we assumed = K and rie = 10^ cm~^, 
appropriate for the values that we have computed 
in the rim (Section 4.2). Theoretical Balmer line 
ratios from Osterbrock (1989) and the reddening 
law of Cardelli, Clayton, & Mathis (1989) were 
used. Our derived value is in fair agreement with 
the average value (0.15) reported by other au- 
thors: Ch/3 = 0.08 ± 0.03 (Hyung & Aller 1995); 
CH/3 = 0.24 ± 0.04 (Lame & Pogge 1996); ch/3 = 
0.20 (Liu et al. 1995); ch/3 = 0.10 (Rubin et al. 
2002). 
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Note that wc applied ch/3 = 0.16 for the fluxes 
of all the different zones under analysis, because of 
the highest S/N ratio of the rim, thus avoiding the 
larger errors in the measured Balmcr line ratios of 
the fainter zones. This assumption of a constant 
extinction throughout the nebula is supported by 
the work of Bohigas, Lopez, & Aguilar (1994), who 
reported that the ch/3 for the different structures 
of NGC 7009 are always lower than 0.28, implying 
small corrections to the observed fluxes. 

4.2. Densities and Temperatures 
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Fig. 2. — Electron densities and temperatures de- 
termined with different line ratios as a function of 
the distance from the center of NGC 7009. The 
positions of the various features are marked in the 
uppermost part of the plot (see Figure 1). Symbols 
are plotted slightly displaced in distance in order 
to avoid overlapping. The horizontal lines in both 
plots give the value of vnvh density/temperature 
estimate for the entire nebula: A^e[S ll](NEB) 
= 4000 ± 300 cm-3 Ne[Cl iii](NEB) = 4500 ± 
400 cm-3; Te[0 iii](NEB) = 10 100 ±7 00 K, and 
Te[N II] (NEB) = 10 300 ± 800 K. 

As described above, eight different regions were 
selected in the spectrum of NGC 7009 (see also 
Figure 1). The extinction-corrected fluxes were 
then used in order to obtain the densities and tem- 
peratures for each of the regions that are plotted 



in Figure 2. These physical parameters were ob- 
tained using the nebular package of IRAF (Shaw 
& Dufour 1995), based on the five-level atom pro- 
gram (De Robertis, Dufour, & Himt 1987) to de- 
rive the physical conditions in a low-density neb- 
ular gas, including diagnostics for a large number 
of ions and emission lines. 

To estimate densities, we assumed an electron 
temperature of 10^ K (see below). In most posi- 
tions, we were able to measure the densities from 
both the [S ii] 6716 A/6731 A and [CI iii] 5517 
A/5537 A line ratios. These provide a simultane- 
ous estimate of densities in low- ionization ([S ll]) 
and higher ionization ([CI III]) regions. In gen- 
eral, we found a good agreement between the re- 
sults from these two density estimators for most of 
the structures (see Table 1 and Figure 2). Using 
the total flux integrated along the slit, we obtain 
iVe[Cl iii]/iVe[S ii] = 1.13; the horizontal dashed- 
dashed line in Figure 2 indicates the A^e[Cl III] 
overall value, while the continuous line shows the 
corresponding Ne[S li] overall density. 

Figure 2 clearly shows that the densities of 
the outermost knots are very similar to those of 
the jets connecting them to the edge of the in- 
ner rim, with A^e(Ki) = 1900 ± 500 cm~^ and 
iVe(Ji)= 1300±200 cm^'' and, on the opposite 
side, A^e(K4) = 1300 ± 200 cm-^ and A^e(J2) = 
1300 ± 500 cm~^. On the other hand, the in- 
ner pair of knots have higher electron densities, 
namely Afe(K2) = 4500 ± 400 cm'^ and A^e(K3) 
= 5000 ± 500 cm~"^, and similar to those of the 
rim (iVe(Ri)= 5500 ± 600 cm-^ and Ne{R2) = 
5900 ± 800 cm-3). 

As for the densities, the temperature estima- 
tors used are appropriate for zones of low- and 
high-excitation. We used the following line ra- 
tios: 7(4959 A + 5007 A)//(4363 A) for Te[0 iilj; 
1(6548 A + 6583 A)//(5755 A) for re[N ii]; and 
/(6716 A/6731 A)/7(4069 A + 4076 A) for Te[S ii]. 
The general trend of temperatures (Figure 2. bot- 
tom) is that, within the errors, they are constant 
throughout the nebula, having an average value of 
re[0 in] = 10 200 K and Te[N ii] = 11 100 K. In 
spite of a small discrepancy of Te [S ii] at the posi- 
tion of R2, the temperature is remarkably constant 
across this nebula. 

There are a number of papers that discuss the 
physical properties of NGC 7009 (Rubin et al. 
2002; Hyung & AUer 1995; Bohigas et al. 1994; 
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Balick ct al. 1994). In general, there is good over- 
all agreement between our temperature determi- 
nations and previous ones. A more detailed zone- 
by-zonc comparison can be made with some of 
these papers. Balick et al. (1994) studied the west 
rim, cap, and ansae. They found Te[N ii] and Te[0 
III] of 10 000 K and 9400 K, respectively, for the 
rim, 9600 K and 9600 K for the cap, and 8100 K 
and 11 500 K for the ansa. With the exception of 
Te[N ll] of that ansa, those figures are very similar 
to ours (see Table 1). From the six regions of the 
PN studied by Bohigas et al. (1994), two coincide 
with ours (those of the W and E ansae), and they 
found, for both structures, an [O III] temperature 
of 9800 K, again consistent, within the errors, with 
our measurements. The Te[0 ill] determined from 
the HST/WFPC2 data by Rubin et al. (2002) is 
also consistent with the values we derived, since 
they found 9000 K < Te[0 ill] < 11000 K. Their 
STIS long-slit results for temperature are limited 
to the 13 central arcseconds along the major axis. 
The average Te[0 ill] for K2, Ri, R2, and K3 are 
~10 300 K ~10 700 K, ~10 600 K, and ~10 100 K, 
respectively, which are very similar to our figures 
in Table 1. 

The above papers also report density estimates, 
whose agreement with our measurements is not, 
however, as good as for the electron temperatures. 
For instance, Balick et al. (1994) found Ne[S 11] 
= 4900 cm-3, 4100 cm'^, and 1000 cm'^ for the 
western rim, cap, and ansae, respectively. They 
also found a large difference between the A''e[S 
ll]/A^e[Cl III] ratio at the position of the W rim 
and cap, which was measured to be ~3 and ^0.4, 
respectively, while ours is ^1.0 and ^0.8 (in line 
with Ne[S II] « jVe[Cl III] fomid by Copetti & 
Writzl 2002 for a large sample of PNe). Lame & 
Pogge (1996) found that the rim plus inner knots 
of NGC 7009 have densities as high as 8000 cm'^, 
and that the rim without including K2 and K3 has 
a significantly lower density, ~4000 cm~'^. There 
is also only moderate agreement between our den- 
sities and those of Bohigas et al. (1994), who ob- 
tained 2300 ± 600 cm-3 and 4300 ± 2700 cm'^ 
for the W and E ansae, respectively. The reason of 
the discrepancies between different authors might 
be the S/N ratio of the relevant lines, and the see- 
ing. Note that our observations are, in general, 
much deeper and were obtained during better see- 
ing conditions than the above studies, allowing a 



more precise definition of the spatial regions. 
4.3. Excitation Properties 
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Fig. 3. — Emission line spatial distribution along 
P. A. = 79°. From top to bottom: the [N Il]/Ha 
(continuous) and [S Il]/Ha (dashed) line ratios; 
the [O ii]/[0 III] line ratio; the [N il] and [O li] 
line profiles for the inner (K2, K3) and outer (Ki, 
K4) pairs of knots. The position of all the selected 
features of Figure 1 are marked at the upper part 
of the plot. 

Concerning the general excitation structure of 
NGC 7009, the low-ionization ([O i], [N 11], [S 11], 
[O 11]) line emission has local maxima at the po- 
sition of the pairs of knots (Ki„4). The higher 
([O III], [Ne III]) ionization emission lines, how- 
ever, are much more prominent within the rim of 
the nebula. 

As the line ratio profiles might provide addi- 
tional interesting information about the excitation 
conditions of a nebula, we show in Figure 3 the [N 
II] /Ha, [S II] /Ha, [O 11] /[O ill], and [N 11], [O ll] 
spatial profiles along the major axis of the nebula. 
The top panel clearly shows the inner and outer 
pairs of knots as prominent zones of low ioniza- 
tion. Moreover, the fact that [N 11] /Ha, [S 11] /Ha 
and [O 11] /[O ill] ratios are very similar in shape 
suggests that their local maxima might not be due 
to the overabundance of a given element. These 
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line ratios arc rather low at the position of the jets 
and the rim implying that the excitation of the jets 
is quite different from that of the knots. The av- 
erage values of [N ii]/Hq! in the outer knots are 
1.00 (Kl) and 0.96 (K4), much higher than those 
usually found for spherical and elliptical PNe (e.g., 
Aller & Liller 1968). The inner pair of knots, on 
the other hand, have less extreme [N 11] /Ha line 
ratios, namely 0.13 and 0.14 for K2 and K3, re- 
spectively. 

The lower two panels of Figure 3 show the [N 
11] and [O 11] spatial profiles at the position of 
the inner (top) and outer (bottom) pairs of knots. 
The [N 11] and [O 11] profiles have been normal- 
ized in the nebular regions just ahead of the knots 
in order to see whether the variation of the [N 
11] and [O 11] emission through the knots shows 
any effects related to coUisional quenching of the 
emission lines. The critical density for coUisional 
de-excitation is 16 000 cm~^ and 3100 cm""^ for 
the [O 11] 3726 A 3729 A doublet, respectively, 
while it is 86 000 cm"-'' for the [N 11] 6583 A line. 
Figure 3 shows that the [N 11] and [O 11] profiles 
are very similar in the outer, lower density (Ng 
< 2000 cm^'') knots, whereas in both inner knots 
(Ne « 5000 cm~^) on the contrary, the [O ll] pro- 
file is much flatter than the [N 11] one. This rep- 
resents a direct proof that coUisional quenching 
is indeed affecting the [O 11] doublet in regions 
of relatively high densities, and caution should be 
exercised when, for example, computing densities 
using the [O 11] doublet. 

Peculiar line ratios from structures like the 
knots of NGC 7009 have been attributed to 
anomalous abundances and/or shock excitation 
(BaUck et al. 1987; Miranda & SoU 1992). In or- 
der to test whether shock excitation plays a role 
in NGC 7009, we show in Figure 4 two relevant 
diagnostic diagrams (c.f. Canto 1981; Phillips & 
Cuesta 1999). The behavior of the different com- 
ponents of NGC 7009 in the Ha/[N 11] (6548 + 
6583) versus Hq/[S ii] (6716 + 6731) diagram is 
very interesting, revealing the eight regions to be 
neatly separated and ordered in pairs along the 
PN band. This is in line with the trends found 
in other PNe (see the case of M 2-9, Phillips & 
Cuesta 1999): the rim appears on the right of the 
diagram, whereas the outer zones of the nebula 
lie in the bottom left of the diagram. An interest- 
ing result from this diagram is that the jets and 
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Fig. 4. — Diagnostic diagrams showing the loci 
of the diff'erent observed structures in NGC 7009: 
Top, log Ra/[N 11] (6548 + 6583) versus log 
Ha/[S II] (6716 + 6731), from Canto (1981). Bot- 
tom, log [O III] (4959 + 5007)/Ha versus log [S 
II] (6716-1-6731) /Ha from Phillips & Cuesta (1999). 



rim, which are high-excitation zones are on the 

right side, while both pairs of knots appear to be 
separated from the former in terms of degree of 
excitation. 

The bottom panel of Figure 4, showing [O 
III] (4959 + 5007)/Ha versus [S ii](6716 + 6731)/Ha, 
separates the zone of radiatively excited emis- 
sion lines (the PN zone, adapted from Phillips 
& Cuesta 1999) from the zone mainly excited by 
shocks (from the predictions of the plane-parallel 
and bow-shock models of Hartigan, Raymond, & 
Hartmann 1987). It is clear that all eight struc- 
tures in NGC 7009 are mainly radiatively excited 
by the PN central star. This is remarkable, since 
other PNe show, at least in some of their mor- 
phological components, indications of shock exci- 
tation, whereas in NGC 7009 not even the jet and 
outer knots approach the shock area in Figure 4. 
On the other hand, the fact that forbidden lines 
are radiatively excited in this nebula validates our 
calculations of physico-chemical parameters for 
the different regions. 
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4.4. Chemical Abundances 
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Fig. 5. — Total abundances profiles for He, O, 
N, Ne and S. The bottom-right panel shows the 
N+/0+ profile. The position of all the features 
selected from Figure 1 are marked in the upper 
part of the plot. 

The chemical composition of a PN results from 
the mixing of the elements produced by the cen- 
tral star and dredged up to the surface with those 
present in the original gas from which the star was 
made. The abundances of O, Ne, Ar, and S are 
thought to represent the chemical composition of 
the interstellar medium when the progenitor star 
was born, since these elements are neither pro- 
duced by the PN progenitor nor very much af- 
fected by nucleosynthesis (see Iben & Renzini 1983 
and Stasihska 2002 for reviews). On the contrary, 
He, C, and N changes dramatically in the progeni- 
tor star enriching the stellar surface and, owing to 
the drcdgc-up episodes, enrich the PN envelopes. 

Abundances of ionized nebulae are calculated 
based either on empirical methods or on model 
fitting (Stasinska 2002). In the former case, com- 
plete multiwavelength data would give the emis- 
sion for all ions of a given element, and the total 
abundance is the sum of the ionic ones. Usually, 
however, only a few lines of some ions can be mea- 
sured and one must correct for unseen ions using 
ionization correction factors. 



Here, we have computed the ionic and total 
abundances (relative to hydrogen) for the different 
regions of NGC 7009, and also for the whole neb- 
ula via the ionization correction factors [icf) fol- 
lowing the scheme in Kingsburgh & Barlow (1994), 
as done, for example, by Corradi et al. (1997). 
The uncertainties in this method have been dis- 
cussed by Alexander & Balick (1997) and appear 
to be significant for the analysis of spatially re- 
solved long-slit observations. 

Table 3 lists the ionic and total abundances and 
the ic/used. Total abundance profiles along the 
major axis of the nebula are shown in Figure 5. 
Errors in the ionic abundances are derived by tak- 
ing into account both the errors in the line ratios 
and those in the adopted temperatures. Errors in 
the total abundances are obtained by propagating 
the errors in the ionic abundances as well as on 
the icj. Typical resulting errors (lo") in the total 
abundances are less than 10% for He (for which 
no icf is applied), 10-20% for O, less than 35% for 
Ne, and 25-40% for N and S. 

First of all note in Figure 5 that the abundances 
obtained for the regions Ji and J2 should be inter- 
preted with great caution because these are indeed 
very faint regions in the nebula. In particular, the 
[N 11] A5755 A fine could not be measured at the 
jet positions, and we could not determine re[N 
11]. Instead, we have assumed re[N 11] = Te[0 ill] 
for Ji and J2, based on the fact that the Te[N il] 
to Te[0 III] ratio of all the other regions ranges 
between 1.01 and 1.25. Therefore, although the 
adopted temperature for Ji and J2 seems reason- 
able, the ionic abundances of the low-ionization 
ions for these regions are probably more uncertain 
than the quoted errors. 

Figure 5 shows that the abundance of Neon re- 
mains constant throughout the nebula. Oxygen 
and Sulfur abundances appear larger in Ki and 
K2 than in the other regions, but the variations 
are under the 2a level. The constancy of these 
elements is expected from stellar evolution theory 
and is consistent with previous results for most 
spatially resolved PNe studied so far (cf. Perinotto 
& Corradi 1998). Helium also remains practically 
constant, implying no or little enrichment in the 
nebular cjecta. 

Figure 5 also shows evidence for an increase of 
the N abundance in the outer knots Ki and K4, as 
previously suggested by Balick et al. (1994). They 
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reported a much higher overabundance, up to a 
factor of 5 for both K4 and K3 (their western ansa 
and cap), although they also warned that their 
uncertainties in the abundances could be as high 
as a factor of 2. 

Our results can also be compared with those of 
other authors, who found "peculiar" abundances 
of N and O for certain zones of NGC 7009. An 
overabundance of N/H in the outer knots of a fac- 
tor of 2 or more with respect to the other regions 
of the nebula was first derived by Czyzak & AUer 
(1979). Baker (1983), on the other hand, found 
that the N abundance does not vary substantially 
through the PN, but that 0/H is marginally (by 
a factor of ~ 2) lower at the position of the outer 
knots. Note, however, that based on much better 
data and an analysis of four PNe containing ansae 
and including NGC 7009, BaUck et al. (1994) con- 
cluded that only N/H increases at the ansae, by 
a factor between 2 and 5 as compared with the 
rim. Such behavior was interpreted as the result 
of the recent ejection of N-enriched high-velocity 
material by the PN central star. 

A good test to investigate further this point is 
the inspection of the profile of the N+ZO"*" ionic 
abundance ratio (commonly taken as a measure of 
the total N/0 abundance that does not depend on 
ic/ corrections). This is shown in the bottom right 
panel of Figure 5, and shows the same trend as the 
N/H ratio. N/0 ratios deduced by Czyzak & AUer 
(1979) also show an apparent gradient from the 
rim to the ansae. Note, however, that Perinotto 
& Corradi (1998) and Hajian et al. (1997) did not 
find any evidence of peculiar N or O abundances in 
the 17 PNe studied, some of which possess ansae. 

Finally, our total abundances of He, O, N, Ne, 
and S for the whole nebula (integrated over the 
entire slit, i.e. those labeled NEB in Table 3) 
are in a general agreement with those obtained by 
Baker (1983), Liu et al. (1995), and Hyung & Aller 
(1995). More accurate abundances were derived 
by Pottasch (2000) combining lUE and ISO/IRAS 
line fluxes with optical data from the literature, 
and are also listed in Table 3. These abundances 
have icf of order of unity and it is therefore rele- 
vant to compare them with our values. Percentual 
differences (from 4% for He to 30-40% for N and 
Ne) are similar to our estimated errors and indi- 
cate a fair level of agreement between the two sets. 

We find that the sulfur abundance for the whole 



nebula is lower than the corresponding solar value, 
while the He, O and Ne ones are approximately 
solar (Grevesse & Anders 1989; Allende-Prieto et 
al. 2001); the N abimdance appears similar or 
marginally higher than the solar value. 

In conclusion, we find a homogeneous elemen- 
tal abundance across the nebula, to within 9%, 
17%, and 35% for He, O, Ne and S, respectively. 
Nitrogen seems to be enhanced in the outermost 
knots of NGC 7009 by a factor < 2, but this evi- 
dence is only marginal considering the errors and 
the enormous range in the derived icf. The un- 
certainties intrinsic to the method are also rather 
large (Alexander & Balick 1997), but the present 
data seem to discard variations at the level found 
by Balick et al. (1994). 

5. Discussions 

How the observed physico-chemical character- 
istics of the microstructures in NGC 7009 fit into 
current models for LIS formation and evolution? 
Let us first note that discrete knots such as Ki 
to K4 could form during the evolution of continu- 
ous jets: dense knots develop at the tips of jets as 
matter accumulates at the head of the jet, while 
knots in the main body of the jet can originate 
via several forms of instabilities (kink, Vishniac, 
Rayleigh-Taylor, and Kelvin-Helmholtz). Some- 
times, both jets and knots are observed, whereas 
in other instances we observe only the latter. From 
a theoretical point of view, the cooling conditions 
on the post-shock gas and the stage of the jet's 
evolution mainly determine whether or not both 
jets and knots are seen (see, for instance, Soker 
1990; Frank, Bahck & Livio 1996 and Garcia- 
Segura 1997). 

According to hydrodynamical (HD) and mag- 
netohydrodynamical (MHD) studies (Garcia- 
Segura et al. 1999: Garcfa-Segura & Lopez 2000; 
Steffen, Lopez, & Lim 2001; Gardiner & Frank 
2001; Blackman et al. 2001), the interplay between 
the stellar AGB and post-AGB winds (for single- 
star progenitors) or between stellar and disk winds 
(if the progenitor is a binary star) can form highly 
coUimated jets in PNe. Jets that originate in this 
way are supersonic, highly coUimated, two-sided, 
and roughly coeval with the main shell (single 
star) or yoimger (binaries). Another character- 
istic of some of the models is the linear increase 
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of the jets' expansion velocity with distance from 
the center (Garcia-Segura et aL 1999; Steffen et 
al. 2001). The orientation of jets that result from 
these models depends on the alignment between 
the shell axis and the disk/rotation (or magnetic) 
axis (Garcia-Segura & Lopez 2000; Blackman et 
al. 2001), as well as on the presence of precession 
or wobbling (Cliffe et al. 1995; Garcfa-Segura 
1997; Livio & Pringle 1997). 

If the kinematical results of Reay & Atherton 
(1985) are correct (but see the discussion in Sec- 
tion 3), then, as we have shown in a previous paper 
(Gongalves et al. 2001), the morphology and kine- 
matics of the jet and knots of NCG 7009 can be 
qualitatively explained in the light of HD or MHD 
single-star interacting stellar wind model. How- 
ever, as far as the comparison of our results with 
models in terms of physical parameters is con- 
cerned, the situation becomes worse. On the one 
hand, most of the above models do not show clear 
predictions for the densities and temperatures of 
jets and knots. On the other hand, those which 
do show such predictions are hardly in agreement 
with the physical parameters derived from the 
present observations, as we will show in the fol- 
lowing. 

The simple process of mass accumulation im- 
plies that knots at the jets' tips are necessarily 
denser than the jets themselves. As we showed 
in Section 4 this is not the case in NGC 7009. 
An additional, very intriguing, property of the 
jets/knots of this PN is that we find no evidence 
for shock excitation from the optical lines (the line 
ratios of the jets and outer knots appear to be es- 
sentially radiativcly excited, see Section 4.3). 

Further insights can be gained by comparing 
our results with the following detailed models. 
Steffen et al. (2001) developed HD stagnation 
knot/jet 2-D models which are based on the slow 
AGB and fast post-AGB winds interaction with 
lack of momentum in the polar direction of the 
main shell. They show that this interaction gives 
rise to "microstructures" which first produce a 
slightly elongated knot attached to the rim, then 
a knot further elongated and detached to the rim, 
and finally a polar jet with a characteristic out- 
ward linear increase in velocity. The stagnation 
knot /jet model has transient lower density jets 
with higher density knots at the tips, but the op- 
posite occurs with the more permanent jet struc- 



tures. Therefore, the stagnation knot /jet mod- 
els do not explain the outer knots and jets of 
NGC 7009, whose outer pair of knots and jets have 
very similar densities, contrary to the model pre- 
dictions. 

Recent 3D MHD models by Garcia-Segura & 
Lopez (2000) can also be qualitatively compared 
with our results. These models seem to be able 
to reproduce successfully the formation of large- 
and small-scale structures in PNc. These au- 
thors state that the mass-loss rate during the 
post-AGB is a key parameter in distinguishing 
betwcion the formation of PNe with highly colli- 
mated jets and PNe with knots (ansae-like struc- 
tures): M > 10"'' Mq yv~^ leads to the forma- 
tion of jets, while lower mass-loss rates (M ~ 
10~^ Mq yr~^) are only able to produce knots, in- 
stead of jets. Although the low resolution of these 
models precludes making quantitative predictions 
for the physical parameters of the knots and jets, 
their results suggest (their Figure 7) that jets and 
knots are indeed shock-dominated while our anal- 
ysis indicates that the jets and the outer knots of 
NGC 7009 are mainly radiatively excited. 

An alternative possibility for the formation of 
polar knots/jets in elliptical PNe is the model 
proposed by Frank et al. (1996). They suggest 
that jets and ansae can be formed during the 
"momentum-conserving" proto-PN stage, i.e., be- 
fore the beginning of the fast wind. The key in- 
gredients of this model are the gradual change 
from a slow to a fast wind, as well as an appre- 
ciable equator-to-pole density contrast in the slow 
wind. The accelerated wind would flow preferen- 
tially parallel to the shock, and the gas would fol- 
low the shocks' prolate boundary toward the polar 
axis. The final result is that the gas flows con- 
verge at the polar axis where they collide, form a 
new shock that redirects the flow into a jet, and 
moves outward along the main axis of the bubble. 
Jets/ansae formed this way can be accelerated up 
to 70kms~^ and have densities > 10^ cm^^, be- 
fore being ionized by the UV photons of the post- 
AGB star. 

In fact the jet/ansae of NGC 7009 look like a 
continuity of its elliptical rim (remembering that 
rim and jet are very similar in excitation degree), 
which is very much in the line of the "converg- 
ing polar flows" model for the formation of jets 
and ansae. However, such flows do not tend 
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to converge into stable "microstructurcs" , since 
the fast/slow wind interaction could destroy these 
structures (Prank et al. 1996; Dwarkadas & Balick 
1998; Lopez 2002). 

Finally, if the nitrogen overabundance of the 
outer pair of knots is taken as real (but see the dis- 
cussion in Section 4.4), the formation of Ki and K4 
might be related to recent high-velocity ejections 
of enriched material from the central star. Such a 
conclusion requires that the high-velocity nature 
of the structures is effectively confirmed by more 
robust kinematical data (especially concerning the 
proper motions). 

Another issue to be further investigated is the 
relation between the thin jet and the "fat" polar 
bubbles of the nebular shell. The jet seems to be 
located along the polar axis of the bubbles, and 
the outer knots lie exactly at the edge of the bub- 
bles. It is then conceivable that the outer knots 
are the results of the interaction of the jet with the 
walls of the bubbles. However, would such interac- 
tion enhance the low-ionization emission without 
developing a significant density contrast, as ob- 
served? 

The formation of the inner system of knots is 
also uncertain. Morphologically, the knots seem 
to be located in the shell, and might well share its 
expansion velocity, as found in other PNe, such as 
IC 2553 and NGC 5882 (Corradi ct al. 2000). If 
so, they could well be the results of in situ instabil- 
ities, or survivors of pre-existing condensations in 
the AGB wind. Our suggestions for the formation 
of the outer and inner pair of knots can readily be 
tested with adequate kinematical modeling of the 
different structures in NGC 7009. 

6. Conclusions 

A physico-chemical investigation of the sym- 
metrical pairs of structures — rim, inner and outer 
knots, and jetlike systems connecting the latter - 
in NGC 7009 is presented. This study, the most 
extensive (eight individual regions along the ma- 
jor axis), better resolved (1 arcsec spatial resolu- 
tion within the central 1 arcmin) and complete to 
date (including all diagnosis from available optical 
lines), has revealed a clear and coherent pattern of 
physical, excitation, and chemical characteristics 
through the nebula: 1) the electron temperature, 
both for low- and high- excitation species, is re- 



markably constant; 2) electron densities are sim- 
ilar for both the jetlike structures and the outer 
knots, but they are roughly one- third of those 
found in the pair of inner knots and the rim; 3) no 
sign of shock excitation is found for any of the mi- 
crostructurcs or the rim; and 4) a notable chemical 
uniformity is found for the eight regions studied, 
with the exception of a marginal, moderate ( x 2) 
enhancement of nitrogen at the outernmost knots. 

The origin of the system of microstructures in 
NGC 7009 is discussed in the light of its physico- 
chemical parameters and the available theoretical 
models, but none of these appears able to account 
for both the lack of a strong density contrast be- 
tween the outer knots and jetlike structures, and 
the absence of shocks in them. More than 200 
years after its discovery by Herschel, the Saturn 
Nebula is still challenging our comprehension. 

7. Acknowledgments 

We would like to thank the referee, Letizia 
Stanghellini, for several useful comments and sug- 
gestions that had help us to improve the paper. 
We are grateful to F. Sabbadin for detecting a 
misidentification (now corrected) of features in 
Pig. 1 of the preprint. The work of DRG, RLMC 
and AM is partially supported by a grant from the 
Spanish Ministry of Science and Technology (AYA 
2001-1646). 

REFERENCES 

Acker, A., Ochsenbein, F., Stenholm, B., Tylenda, 
R., Marcout, J., & Schohn, C, 1992, Strasbourg- 
ESO catalogue of Galactic planetary nebulae, 
ESQ 

Alexander, J., & Balick, B. 1997, AJ, 114, 713 

AUende-Prieto, C, Lambert, D. L., & Asplund, 
M. 2001, ApJ, 556, 63 

Aller, L. H., & Liller, W. 1968, in Nebulae and 
Interstellar Matter, ed. B. M. Middlehurst & 
Aller (Chicago: Univ. Press), 498 

Baker, T. 1983, ApJ, 267, 630 

Balick, B. 1987, AJ, 94, 671 

Balick, B., Alexander, J., Hajian, A. R., Terzian, 
Y., Perinotto, M., & Patriarchi, P. 1998, AJ, 
116, 360 



11 



Balick, B., Owen, R., Bignell, C. R., & Hjellming, 

R. M. 1987, AJ, 94, 948 

Balick, B., Perinotto, M., Maccioni, A., Alexan- 
der, Terzian, Y., & Hajian, A. R. 1994, ApJ, 

424, 800 

Balick B., Preston, H. L, & Icke, V. 1987, AJ, 94, 
164 

Balick B., Rugers, M., Terzian, Y., & Chengalur, 

J. N. 1993, ApJ, 411, 778 

Blackman, E. G., Frank, A., & Welch, C. 2001, 
ApJ, 546, 288 

Bohigas, J., Lopez, J. A., & Aguilar, L. 1994, 
A&A, 291, 595 

Calm, J. H., Kaler, J. B., & Stanghellini, L. 1992, 

A&AS, 94, 399 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 
1989, ApJ, 345, 245 

Canto, J. 1981, In: Investigating the universe, by 
D. Reidel Publishing Co., p. 95 

Cliffe, J. A., Frank, A., & Livio, M., Jones, T. W. 
1995, ApJ, 447, L49 

Copetti, M. V. F., & Writzl, B. C. 2002, A&A, 

382, 282 

Corradi, R. L. M., Manso, R., Mampaso, A., & 
Schwarz, H. E. 1996, A&A, 313, 913 

Corradi, R. L. M., Gongalvcs, D. R., Villaver, E., 
Mampaso, A., & Perinotto, M. 2000, ApJ, 542, 

861 

Corradi, R. L. M, Perinotto, M., Schwarz, H. E., 
Claeskens, J-F 1997, A&A, 322, 975 

Czyzak, S. J., & Aller, L. H. 1979, MNRAS, 188, 
229 

Dwarkadas, V. V., & Balick, B. 1998, ApJ, 497, 
267 

De Robertis, M., Dufour, R. & Hunt, R. 1987, 
JRASC, 81, 195 

Frank, A., Balick, B., & Livio, M. 1996, ApJ, 471, 
L53 

Frank, A., & Mellema, G. 1994, ApJ, 430, 800 

Garcfa-Segura, G. 1997, ApJ, 489, 189 

Carci'a-Sogura, G., Langcr, N., Rozyczka, M., & 
Franco, J. 1999, ApJ, 517, 767 

Garcia-Segura, G., & Lopez, J. A. 2000, ApJ, 544, 
336 



Gardiner, T. A., & Frank, A. 2001, ApJ, 557, 250 

Goncalvcs, D. R., Corradi, R. L. M., & Mampaso, 
A. 2001, ApJ, 547, 302 

Gorny, S. K., Schwarz, H. E., Corradi, R. L. M., 
& van Winckel, H. 1999, A&AS, 136, 145 

Grevesse, N., & andcrs, E. 1989, in Waddington J., 
ed.. Cosmic Abundances in Matter, AIP, New 
York, p. 1 

Guerrero, M. A., Gruendl, R. A., & Chu, Y.-H. 
200 2, A&A, 387, LI 

Hajian, A. R., Balick, B., Terzian, Y., & Perinotto, 
M. 1997, ApJ, 487, 313 

Hartigan, P., Raymond, J., & Hartmann, L. 1987, 
ApJ, 316, 323 

Hyung, S., & Aller, L. H. 1995, MNRAS, 273, 958 

Ibcn, I, & Rcnzini, A. 1983, ARA&A, 21, 271 

Kingsburgh, R. L., & Barlow, M. J. 1992, MN- 
RAS, 257, 317 

Kingsburgh, R. L., & Barlow, M. J. 1994, MN- 
RAS, 271, 257 

Kwok, S., Purton, C. R., Fitzgerald, P.M. 1978, 
ApJ, 219, L125 

Lame, N. J., & Pogge, R. W. 1996, AJ, 111, 2320 

Liller, W. 1965, PASP, 77, 25 

Liu, X.-W., Storey, P. J., Barlow, M. J., Clegg, R. 
E. S. 1995, MNRAS, 272, 369 

Livio, M. & Pringlc, J. E. 1997, ApJ, 486, 835 

Lopez, J. A., Vazquez, R., & Rodriguez, L. F. 
1995, ApJ, 455, L63 

Lopez, J. A. 1997, in Planetary Nebulae, Proceed- 
ings of the 180th Symposium of the lAU, ed. H. 
J. Habing. and H. J. G. L. Lamers. Publisher: 
Dordrecht Kluwer Academic Publ., p. 197 

Lopez, J. A. 2002, in "Emission Lines from Jet 
Flows", Eds. W.J. Henney, W. Steffen, A.C. 
Raga & L. Binette, RevMexAA (Conference Se- 
ries), 13, p. 139 

Maciel, W. J. 1995, Ap&SS, 229, 203 

Manchado, A., Guerrero, M. A., Stanghellini, L., 
& Serra-Ricart, M. 1996, The lAC Morpholog- 
ical Catalog of Northern Galactic PNe, lAC 
Publ. 

Mellema, G., & Frank, A. 1995, MNRAS, 273, 401 



12 



Mendez, R. H., Kudritzki, R. P., & Herrero, A. 

1992, A&A, 260, 329 

Miranda, L. F., & Solf, J. 1992, A&A, 260, 397 

Osterbrock, D. E., 1989, "Astrophysics of Gaseous 
Nebulae and Active Galctive Nuclei". Univer- 
sity Science Books, cd. Mill Valcy, California. 

Perinotto, M. & Corradi, R. M. L. 1998, A&A, 
332, 721 

Phillips, J. P., & Cueta, L. 1999, AJ, 118, 2919 

Pottasch, S. R. 2000, in Asymmetrical Planetary 
Nebulae II: From Origins to Microstructures, 
ASP Conference Series, Vol. 199., Edited by J. 
H. Kastner, N. Soker, and S. Rappaport, p. 289 

Rcay, N. K., & Atherton, P.D. 1985, MNRAS, 215, 

233 

Rubin, R. H., et al. 2002, MNRAS, 334, 777 

Shaw, R. A., & Dufour, R. J. 1995, PASP, 107, 
896 

Schwarz, H. E., Corradi, R. L. M., & Melnick, J. 

1992, A&AS, 96, 23 

Schonberner, D. 2002, in "Planetary Nebulae and 
Their Role in the Universe" , proceedings of lAU 
Symposium No. 209, ASP Conference Series 

Soker, N. 1990, AJ, 99, 1869 

Stasihska, G. 2002, to be published in the proceed- 
ings of the XIII Canary Islands Winter School 

of Astrophysics (astro-ph/0207500) 

Steffen, W., Lopez, J. A., & Lim, A. 2001, ApJ, 

556, 823 

Zhang, C. Y., & Kwok, S. 1998, ApJS, 117, 341 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v5.0. 



13 



Table 1 

Measured line fluxes normalized to H/3=100. 
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Hcii 4685.7 
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2.62 


10.2 


25.6 


23.0 


8.42 


2.16 


1.21 


15.6 


[Ariv] + Hel 4711+13 


2.16 


4.62 


3.23 


5.40 


4.89 


2.78 


5.08 


1.97 


4.04 


[Ariv] 4740.2 


1.13 


3.63 


2.78 


5.50 


4.49 


2.54 


3.95 


1.48 


3.88 


H/3 4861.3 


100 


100 


100 


100 


100 


100 


100 


100 


100 


[Oiii] 4958.9 


428 


473 


435 


405 


428 


424 


458 


454 


421 


[Oiii] 5006.86 


1256 


1339 


1268 


1179 


1243 


1227 


1381 


1329 


1223 


[Arm] 5191.8 






0.09 












0.07 


[Ni] 5200.2 


6.06 




0.06 










2.82 


0.11 


[Fern] 5270. 


















0.08 



Table 1 — Continued 







J 1 










To 


Ka 


NEB 


[Feiii] + [Cav 5292. 


_ 


_ 


_ 


_ 


_ 


_ 


_ 


_ 


0.09 


Hen 5411.6 


_ 


_ 


0.72 


1.48 


1.67 


0.95 


_ 


_ 


1.12 


[Cliii] 5517.7 


_ 


_ 


0.71 


0.46 


0.46 


0.55 


0.91 


1.03 


0.58 


[Cliii] 5537.9 


_ 


0.40 


0.85 


0.57 


0.59 


0.71 


0.79 


0.96 


0.68 


[Oi] 5577.4 


0.51 




0.14 










. 


0.01 


[Oiii] 5592.4 


. 






0.03 


0.40 








0.05 


[Nil] 5754.6 


6.95 


_ 


0.72 


0.15 


0.20 


0.69 


_ 


4.19 


0.50 


Hei 5875.7 


20.3 


17.0 


17.9 


15.0 


15.3 


15.8 


15.8 


16.2 


15.7 


[Kiv] 6102. 


_ 


_ 


_ 


0.16 


0.18 


0.11 


_ 


_ 


0.15 


Hen 6166 + 6170 


















0.06 


Heii 6234. 


















0.06 


[Oi] 6300.3 


31.0 


_ 


1.02 


_ 


_ 


1.69 


_ 


14.8 


1.06 


[Siii] 6312.1 


4.29 


1.80 


2.84 


1.41 


1.26 


1.97 


1.50 


3.59 


1.86 


Sill 6347.1 


_ 


_ 


_ 


0.07 


_ 


_ 


_ 


_ 


_ 


[Oi] 6363.8 


11.3 


_ 


0.38 


0.12 


0.17 


0.56 


_ 


5.33 


0.36 


Sin 6371. 


_ 


_ 


_ 


0.10 


0.12 


_ 


_ 


_ 


_ 


Heii 6406.5 








0.19 


0.19 










[Nil] 6548.0 


132 


10.8 


17.9 


3.57 


3.57 


13.9 


7.05 


72.5 


11.2 


Ha 6562.8 


394 


337 


408 


312 


341 


280 


217 


224 


313 


[Nil] 6583.4 


397 


19.5 


53.8 


7.93 


6.26 


39.2 


14.7 


217 


30.2 


Hei 6678.1 


8.41 


5.52 


6.16 


4.35 


4.44 


4.00 


3.74 


3.57 


4.46 


[Sii] 6716.5 


41.4 


2.52 


4.82 


0.54 


0.43 


3.13 


2.59 


26.1 


2.63 


[Sii] 6730.8 


57.4 


3.09 


8.19 


0.95 


0.78 


5.42 


3.05 


32.0 


4.34 


Ne[S II] 


2000 


1300 


4500 


5500 


5900 


5000 


1100 


1300 


4000 


iVe[Cl III] 






4700 


5200 


5900 


6000 


1300 


1900 


4500 


Te[0 III] 


9600 


10400 


9300 


10000 


10200 


10100 


11600 


10400 


10100 


Te[N II] 


11000 




9400 


10400 


12 800 


10400 




11700 


10300 


Te[S II] 


7100 




8300 










9400 





Table 2 

Estimated percentage errors in line fluxes in Table 1 and H/3 fluxes 



Line fliixes 


Ki 


Ji 


K2 


Ri 


R2 


K3 


J2 


K4 


NEB 




4.9'^ 


4.2* 


8.4* 


3.5* 


1.4* 


9.1* 


4.9* 


4.9* 


55.3* 


(0.01-0.05)Ih/3 


26 


31 


5 


6 


10 


4.5 


11 


13.5 


4 


(0.05 0.15)Ih/3 


14 


10 


3.5 


4 


6.5 


3.5 


7.5 


10 


3.5 


(0.15-0.30)Ih/3 


9 


8 


3.5 


3.5 


4.5 


3.5 


6 


7.5 


3 


(0.30-2.0)Ih/3 


6 


6 


3.5 


3.5 


3.5 


3 


4.5 


5.5 


3 


(2.0-5.0)Ih;3 


5 


4.5 


3 


3.5 


3.5 


3 


4 


5 


3 


(5.0-10.0)Ih/3 


4.5 


4 


3 


3 


3 


3 


4 


5 


3 


> 10 Ih/3 


4 


4 


3 


3 


3 


3 


4 


5 


3 




1.05(6) 


1.24(6) 


30.2(4) 


30.7(4) 


15.5(4) 


70.2(3) 


1.60(5) 


1.28(6) 


188(3) 



^Sizc (in arcscc) of the selected structure along the slit. 

''In units of 10~'^ erg cm~^ s~'. Values within brackets are the estimated percentage errors in the H/3 
line emission. 



Table 3 
Ionic/total abundances^ 





Ki 


Ji 


K2 


Ri 


R2 


Ks 


J2 


K4 


NEB 


NEB - Pottasch 


He+/H 


1.02E-1(06) 


1.23E-1(05) 


1.04E-1(04) 


9.54E-2(04) 


9.18E-2(04) 


9.46E-2(04) 


1.07E-1(05) 


9.55E-2(06) 


9.77E-2(04) 




He^+/H 


- 


1.2E-3(10) 


5.02E-3(06) 


1.24E-2(07) 


2.43E-2(06) 


1.07E-2(07) 


1.87E-3(08) 


- 


1.30E-2(06) 




He/H 


1.02E-1(09) 


1.24E-1(0S) 


1.09E-1(06) 


1.08E-1(08) 


1.16E-1(06) 


1.05E-1(06) 


1.08E-1(08) 


9.55E-2(09) 


1.11E-1(06) 


1.107E-1 


0°/H 


4.5E-5(11) 


- 


2.40E-6(09) 


5.0E-7(13) 


3.1E-7(16) 


2.61E-6(08) 


- 


1.8E-5(13) 


1.71E-6(08) 




0+/H 


7.5E-5(19) 


2.2E-5(13) 


2.4E-5(14) 


7.0E-6(15) 


1.3E-6(16) 


1.4E-5(14) 


6.06E-6(17) 


4.2E-5(21) 


1.2E-5(13) 




\J /in 








A 1 9 "FT A f <~\'7\ 










/I 91 T? /I 




icf{0) 


1.00 


1.00 


1.03 


1.08 


1.16 


1.07 


1.01 


1.00 


1.08 




O/H 


5.8E-4(15) 


8.4E-4(12) 


6.28E-4(09) 


4.5E-4(11) 


4.82E-4(09) 


4.7E-4(11) 


3.0E-4(15) 


4.5E-4(17) 


4.71E-4(09) 


5.2E-4 


N^/H 


8.7E-6(16) 




2.2E-7(12) 










2.7E-6(19) 


2.7E-7(11) 




N+/H 


5.00E-5(08) 


3.09E-6(07) 


9.66E-6(06) 


1.10E-6(07) 


5.19E-7(08) 


5.59E-6(06) 


1.7E-6(10) 


2.38E-5(09) 


4.45E-6(06) 






7.8 


20.1 


25.7 


64.7 


355 


32.1 


50.5 


10.6 


38.2 




N/H 


3.8E-4(35) 


6.2E-5(26) 


2.4E-4(24) 


r.0E-5(26) 


1.8E-4(29) 


1.8E-4(24) 


8.6E-5(34) 


2.5E-4(38) 


l.rE-4(23) 


1.3E-4 


Ne^+ /H 


9.8E-5(14) 


1.1E-4(12) 


1.1E-4(10) 


1.0E-4(10) 


9.6E-5(10) 


1.0E-4(10) 


9.3E-5(17) 


1.2E-4(15) 


1.05E-4(09) 




tcf(Ne) 


1.14 


1.06 


1.07 


1.10 


1.17 


1.11 


1.03 


1.10 


1.11 




Ne/H 


l.lE-4(29) 


l.lE-4(26) 


1.2E-4(20) 


1.1E-4(21) 


l.lE-4(20) 


l.lE-4(20) 


9.6E-5(34) 


1.3E-4(30) 


1.1E-4(18) 


1.9E-4 


S+/H 


2.2E-6(13) 


8.!)9E-8(08) 


5.11E-7(08) 


4.4E-8(11) 


2.2E-8(13) 


2.7E-7(10) 


7.0E-8(10) 


1.0E-6(14) 


2.09E-7(06) 




S^+/H 


7.45E-6(21) 


3.2E-6(15) 


7.4E-6(14) 


2.1E-6(16) 


9.7E-7(19) 


3.4E-6(14) 


1.6E-6(19) 


4.8E-6(23) 


3.3E-6(13) 




icAS) 


1.43 


1.91 


2.07 


2.79 


4.91 


2.22 


2.58 


l..';7 


2. 3. 1 




S/H 


1.39E-5(30) 


6.3E-6(23) 


1.6E-5(23) 


6.1E-6(26) 


4.9E-6(30) 


8.1E-6(23) 


4.4 E-6(29) 


9.3E-6(35) 


8.3E-6(26) 


ll.lE-6(26) 



Estimated percentage errors in brackets. 



